High-energy neutrino conversion and the lepton asymmetry in the universe 



CLunardini ), A.Yu.Smirnov 6 ) 
a) SISSA-ISAS, via Beirut 2-4, 34100 Trieste, Italy 
and INFN, sezione di Trieste, via Valerio 2, 34121 Trieste, Italy 

b) The Abdus Salam ICTP, Strada Costiera 11, 34100 Trieste, Italy 
and Institute for Nuclear Research, RAS, Moscow, Russia 

We study matter effects on oscilfations of high-energy neutrinos in the Universe. Substantiai 
effect can be produced by scattering of the neutrinos from cosmoiogical sources (z Ji 1) on the relic 
neutrino background, provided that the latter has large CP-asymmetry: r\ = (n v — n&)/n 7 <; 1, 
where n„, and n 7 are the concentrations of neutrinos, antineutrinos and photons. We consider in 
details the dynamics of conversion in the expanding neutrino background. Applications are given to 
the diffuse fluxes of neutrinos from GRBs, AGN, and the decay of super-heavy relics. We find that 
the vacuum oscillation probability can be modified by ~ (10 — 20)% and in extreme cases allowed by 
present bounds on r\ the effect can reach ~ 100%. Signatures of matter effects would consist (i) for 
both active-active and active-sterile conversion, in a deviation of the numbers of events produced 
in a detector by neutrinos of different flavours, N a (a = e, fi, r), and of their ratios from the values 
given by vacuum oscillations; such deviations can reach ~ 5 — 15%, (ii) for active-sterile conversion, 
in a characteristic energy dependence of the ratios Ne/N^, N e /N T , iV^/iVV. Searches for these matter 
effects will probe large CP and lepton asymmetries in the universe. 

PACS: 14.60.Pq,13.15+g,98.70.Sa,98.70.Vc 



I. INTRODUCTION 



The detection of high-energy cosmic neutrinos and detailed studies of their properties are among the main challenges 
in astrophysics and cosmology. This will give unique information about the structure of the universe, mechanisms of 
particle acceleration, sources of cosmic rays, properties of the galactic and intergalactic media. They will have also 
important implications for neutrino properties (masses, mixings, etc.) and for particle physics in general. 

Intense fluxes of neutrinos, with energies up to ~ 10 21 eV, are supposed to be produced by cosmoiogical objects 
like Gamma Ray Bursters (GRBs) and Active Galactic Nuclei (AGN) |[] . R was suggested that neutrinos of energies 
as high as 10 22 — 10 24 eV could be produced by topological defects like cosmic strings, necklaces and domain walls 
Furthermore, neutrinos produced by the decay of super-heavy particles have been considered in connection to the 
problem of ultra-high energy cosmic rays exceeding the GZK cutoff j|. 

The detection of high-energy neutrinos from cosmoiogical sources is challenging for neutrino telescopes. The existing 
large water, ice, or airshower experiments open some possibility of detection; however a detailed study requires 
larger detectors to be realized in future. 

The properties of high-energy neutrino fluxes can be modified by oscillations on the way from the sources to 
the Earth. In particular it was marked Q that oscillations lead to the appearance of tau neutrinos in the high- 
energy neutrino flux. Moreover, the study of oscillation effects opens the possibility to probe neutrino mixings and 
distinguish between different mass spectra |(|. In all these studies vacuum oscillations have been considered only. In 
this connection, we address here two questions: 

1. Are matter effects important for high-energy neutrinos propagating in the universe? 

2. Which information on the properties of the interstellar and intergalactic medium can be obtained from the study 
of these effects? 

During their travel from the production point to the detector, the neutrinos cross large amounts of matter, which 
could induce significant refraction and conversion. In ref. we considered the interaction of neutrinos with the matter 
of the source for neutrinos produced in GRBs and AGN. The effects of matter on vacuum oscillations appeared to 
be small. It was also found j?]] that the neutrino-neutrino interaction in the dark matter halos of galaxies does not 
affect the vacuum oscillations significantly Conversely, strong matter effects arc not excluded for neutrinos crossing 
media of larger size, like the halos of clusters of galaxies. Furthermore, neutrinos from cosmoiogical sources travel for 
so large distances in the intergalactic space that the universe itself, with its particle content, can be considered as a 
medium producing refraction effects. In H we found that significant conversion can be realized for neutrinos crossing 
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cosmological distances in the universe with strongly CP-asymmetric neutrino background. 



In this paper we analyze this possibility in detail. We discuss the refraction and conversion effects of the background 
on high-energy neutrinos from cosmological sources and on neutrinos of the background itself. 

Let us describe the relic neutrino gas by the number densities of the various flavours, n a (a — e, e, fi, etc.), and by 
the CP-asymmetry r\ v defined as: 

r) v = (n a - n & )/n~, , (1) 

where n 7 is the concentration of photons. 

The Big Bang Nucleosynthesis (BBN) and structure formation J8|,^| admit large CP-asymmetries for muon and tau 
neutrinos, while the asymmetry for the electron neutrino is strongly constrained: 

|^, r |<10, -0.01 < Ve < 0.3 . (2) 

Large asymmetries have also important implications on the properties of the spectrum of the cosmic microwave 
background radiation (CMBR]_H. The recent results on the second acoustic peak of the CMBR from BOOMERANG 
and MAXIMA-1 experiments |io|| seem to favour a large lepton asymmetry, r\ v ~ 1 |[l]-[l3|. In particular, a satisfactory 
interpretation of the data requires [Q : 

|t? w |£2.2, (3) 

thus providing a stronger restriction of the allowed ?7 M and r\ T values with respect to @) . In our discussion we will 
consider asymmetries 7y M and rj T as large as \r}^^\ ~ 1 — 2, according to the upper limit (^|); however results will be 
given also for larger values, allowed by the less stringent bound (|^). 

We want to underline here that the realization of large CP-asymmetries in the individual lepton flavours is consistent 
with zero lepton asymmetry. This corresponds to: 

Tie + + Vr = , (4) 

that is, to zero total lepton number. Large lepton asymmetry, in contrast, implies large CP-asymmetry. 

The paper is organized as follows. In section |H| the properties of the relic neutrino background are discussed. We 
show that significant matter effects on high-energy neutrinos require large CP-asymmetry of the background, and 
study the refraction and conversion effects in the background itself. In sections [II and |v| we study matter effects on 
high energy neutrinos produced at cosmological distances. Applications are given in section [v| to the diffuse fluxes of 
neutrinos from AGN, GRBs and from the deca y of heavy relics. In section VI we discuss the experimental signatures 



of matter effects. Conclusions follow in section VII 



II. ACROSS THE UNIVERSE 



Let us consider the interactions of high-energy neutrinos propagating from cosmological sources to the Earth. These 
neutrinos cross layers of matter in the source itself, then interact with particles in the interstellar and intergalactic 
media, and finally interact in the matter of our cluster of galaxies and of our galaxy. 

In what follows we will discuss interactions in the intergalactic medium. The effects of the matter of the sources can 
be neglected 0. As we will show later, for neutrino oscillation parameters and energies relevant for this discussion 
also the effect of the galactic halo and of the halo of the cluster of galaxies are very small m . 



A. Minimum width condition 

The necessary condition for significant matter effect is the minimum width condition Q. We define the width of 
the medium as the integrated concentration of the particle background of the universe along the path travelled by 
the neutrino beam: 

d = J ° n{t')dt' . (5) 

Here t is the epoch of production of the neutrinos and to ~ 1 18 s is the present epoch. Using the scaling of the 
concentration n(t) oc t~ 2 , one finds: 
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d(t) = du 
du = t n , 



to 
t 



du 



(l + z)*-l 



(6) 
(7) 



where no is the present concentration of the background and we have introduced the redshift z = (to A) 
The minimum width condition can be written as [ffll: 



1. 



r = — > 1 
do 



where do is the refraction width of the medium: 



do 



7r n 
2V 



(8) 



(9) 



with V being the effective matter potential in a given neutrino conversion channel. The width do corresponds to the 
distance in matter at which the oscillation phase induced by matter equals it/2. Notice that, since V oc n, do does 
not depend on the density of the medium and is determined by the properties of the interaction. The usual weak 
interaction gives do ~ 1/Gf, where Gf is the Fermi constant. 

Let us consider the fulfillment of the condition (|^) for different components of the intergalactic medium. 

1). Due to the very small concentration of nucleons and electrons, the width of these components is extremely small, 
ds/do -C 1, even for neutrinos produced at cosmological distances. Indeed, the baryon concentration can be estimated 
as Hb = rijrjB where t\b = 10 -10 — 10 -9 is the baryon asymmetry of the universe and n 7 the photon concentration. 



At present time n T 
r = r B ~ lO" 11 . 



412 cm 3 . Taking, e.g., production epoch z~lwe find from eqs. that for baryons 



2) . The scattering on the electromagnetic background has negligible effect due to the smallness of interaction. The 
neutrino-photon potential is of the second order in the Fermi constant and depends on the energy of the neutrino 
beam and on the temperature and concentration of the photon gas 14,1^|. Using the results of ref . EH] (see also the 
discussion in p6|) we find from eqs. 

®-® r i < 10_S for neutrino energy E % 10 21 eV and production epoch z ~ 1. 

3) . The effect of the scattering on the neutrino background can produce significant effect if the background has large 
CP-asymmetryQ. Indeed, for asymmetry 77^ ~ 1 we get d u ~ do- Clearly, if the lepton asymmetry is of the order of 
the baryon one, r\ v ~ r)B, the width is negligibly small: d v ~ ds <C do- 
Let us consider the minimum width condition for neutrino background in more detail. The effective potential due 

to the scattering of neutrinos on the relic neutrino background can be written as 



V = Fr)uV2G F r 



'7 } 



(10) 



where F is a constant of order 1 which depends on the specific conversion channel (see sections III A and [V A). With 
the potential ©, eqs. (|), (§) and (§) give the condition: 



1.6- lO -2 !^!^ (z + l)i - 1 



> 1 



For neutrinos produced in the present epoch, z ~ and values of 77^ allowed by the bounds (0)-(| 
is not satisfied: r(z ~ 0) -C 1. From (fTlJ) we can define the epoch Zd which corresponds to r = 1: 



(11 

the condition dl 



Zd = 



1 



1.6 ■ io- 2 !^ 



(12) 



1 For CP-symmetric background significant effects can appear at large temperature, T ^ 1 MeV, due to thermal effects 
or at extremely high neutrino energies, due to neutrino-antineutrino scattering in the resonant Z° channel M. 
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so that for neutrinos produced at z > Zd the minimum width condition is fulfilled. Taking for instance r\ v = 1 and 
F = 2 we get z c i ~ 9; by requiring r ~ 0.3 (which corresponds to 10% matter effect J?J) we find ~ 3.7. 

The following remark is in order. The condition (|^) is necessary but not sufficient to have significant matter effects. 
In particular, for the case of oscillations in uniform medium and small mixing angles (sin 29 0.3) we have found in 
ref. (?j that the width d m i n needed to have conversion probability larger than 1 /2 equals: 

''" >d . (13) 



tan26» 



This quantity represents an absolute minimum. For media with varying density the required width is larger than 

d>min ■ 

We conclude, then, that the only component of the intergalactic medium which can produce a significant matter 
effect is a strongly CP-asymmetric neutrino background, with r\ v ^ 1. Moreover, cosmological epochs of neutrino 
production are required: z 3. 

B. Properties of the relic neutrino background 

Neutrino mixing and oscillations modify the flavour composition of the neutrino background, so that one expects 
the present values of the CP-asymmetries in the various flavours to be different from those at the epoch of BBN; the 
latter are constrained by the bounds (||). 

In this section we assume that large CP-asymmetries are produced at some epoch before the BBN, i.e., at tempera- 
ture T ^ Tbbn — 1 MeV, and study how they evolve with time. The evolution of the flavour densities n e , and n T 
is a non-linear many-body problem, which, in general, requires a numerical treatment |l8| , |l9| ]. In some specific cases, 
however, an analytical description is possible (gCJ and conclusions can be obtained on general grounds. 

1. Three-neutrino system evolution 

Let us first consider the case of mixing between three active neutrinos, v e , v r . In refs. p^j20[ | it has been shown 
that the evolution of the flavour densities has peculiar aspects for the ideal case of a monoenergetic gas of neutrinos 
(with no antincutrinos, n e — 0) initially produced in flavour states. For this specific ensemble of neutrinos the po- 
tential due to neutrino-neutrino interaction cancels in the evolution equation, so that the collective behaviour of the 
system is described by vacuum oscillations. This result holds with a good approximation Jl8| also for realistic neutrino 
energy spectra and in presence of a small component of antineutrinos. For this reason it can be applied to our case of 
interest, in which the background is strongly CP-asymmetric, n„ 3> n e , and neutrinos have a thermal spectrum. In 
what follows we approximate the neutrino energies with the average thermal energy of the gas: E ~< E v >= aT„, 
where T v denotes the temperature of the neutrino gas. The numerical factor a depends on the CP-asymmetry of the 
background: we have a ~ 3.15 in absence of asymmetry, rj ~ 0, and a ~ 3.78 for 77 ~ 1. 

The length scale of flavour conversion is given by the vacuum oscillation length: 

AttE 4na(3T 



Am 2 Am 2 

.5 



2.48 • 10° cm a(3 — — - „ 

H \ lMeV J \ Am 2 



T \ /10" 3 eV 2 



(14) 



where j3 is the ratio between the temperature of the neutrino background and the temperature T of the electromag- 
netic radiation: (3 = T v /T. We have = 1 before the electron-positron recombination epoch, T <; 0.5 MeV, and 
(3 = (4/11) 1 / 3 after this epoch. 

Besides oscillations, for temperatures T <; 1 MeV other phenomena, and therefore other length scales, are relevant: 

• Inelastic collisions. Let us consider a system of two mixed neutrinos, v a , i/j, (a, b = e, v, r). After its production 
as a flavour state, e.g. v = v a , a neutrino oscillates in vacuum until a collision occurs with a particle X of the 
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background. At the time of the collision the quantum state of the neutrino is a coherent mixture of the two 
flavours: v — Av a + Bv^- The effects of the collision depend on the specific reactions that take place |fjl|-fj3f (see 
also the discussion in If the reaction is a scattering, vX — > vX, and the interaction is flavour blind, i.e. 

it is the same for the two flavours a and b, after the collision the neutrino continues to propagate as a coherent 
superposition of v a and v\, and the collision does not affect oscillations. For scattering with flavour-sensitive 
interaction or for absorption processes^, vX —* any, the effect of the collision is to break the coherence between 
v a and Vb, so that after the collision the two flavours evolve independently, developing vacuum oscillations until 
the next collision happens. 

As a result, one easily obtains that for a beam of neutrinos propagating in a medium oscillations are damped 
according to the expression: 



n a {L) 



exp 



L 
X 



In 



1 



1 - 2P r 



(15) 



where n Q a and n a (L) are the fractions of v a in the neutrino beam at the production time and at distance L from 
the production point. Here l c is the coherence length, which represents the distance between two collisions, 
and P c is the vacuum oscillation probability between two collisions: P c ~ sin 2 26 sm 2 (irl c /l v ). From eq. ( |T5| ) 
we see that, if P c ^ 0, with the increase of L (i.e. of the number of collisions) the interplay of oscillations and 
collisions leads to the equilibration of the flavour densities: n a {L — > oo) = n^iL — > oo) = 1/2. The convergence 
to this limit is determined by the equilibration length l eq = l c j In (1/|1 — 2P C |). For small conversion probability, 
P c <C 1, the length l eq is much larger than l c : l eq ~ l c j (2P C ) 3> l c - This is the case if the vacuum mixing is small 
and/or the collisions are much more efficient than oscillations, l c <C l v , so that the vacuum oscillation phase 
$ = 2itl c /l v is small. Thus, equilibration of the flavour densities can be obtained only after a large number of 
collisions: n co a = L/l c 1/(2P C ). Conversely, if P c ~ 1 equilibration is achieved rapidly after few collisions: 
n C oii <^ L/l c . This circumstance is realized if l c ^ l v and the mixing is large, sin 2 29 ~ 1. 

Let us consider the coherence length, l c , in more detail. According to refs. ]2l|,p2|| l c can be written as: 



lc{a,b)~ 



1 



[T abs (a) +T abs (b) +T fs (a,b)] , 



(16) 



where T abs (x) is the rate of absorption processes for the neutrino of flavour x and T^ s (a, b) is the contribution 
of the flavour sensitive scatterings vX — » vX. This quantity is determined by the square of the difference of the 



v a -X and Vy-X scattering amplitudes |2l|2^]. In terms of the total scattering rates T(a) and T(b) one gets pl| : 

Tf s (a, b) ~ r(o) + r(6) - 2 y/T(a)T(b) . (17) 
From cqs. (^6|)-(|i7|), using the rates given in rcf. [24[] we find: 

l c (a, b) = [k(a, b)G 2 F a 2 (3 5 T 5 ] _1 , (18) 



where, for T > 0.5 MeV: 
k(e,n) 



6.5 • 10- 3 [16 + 0.5 (r(Q + r(^)) + 5.3r(-^) + 3.5r(-^)] , 
: 6.5 ■ 10- 3 [0.5(r(f M ) + r(^)) + 3.5(r(-^) + r (-&))] . 



(19) 
(20) 



Here r(£) = I(£)/I(0) and /(£) = J °° x 2 dx/(l + exp(x — ^)). We define £ = n/T with fi the chemical potential 
of the neutrino gas^. At T ~ 0.5 MeV electrons and positrons annihilate, so that for T < 0.5 MeV their 
contribution to the scattering and absorption rates becomes negligibly small and we get: 



2 In the situation we are considering the neutrinos are in thermodynamical equilibrium, so that their disappearance through 
a given reaction is balanced by their production through the inverse process. 

3 Let us recall the relation between the quantity £ and the CP-asymmetry r\: r\= (£ 3 + 7r 2 £)/3 3 /(12£(3)). In eqs. ([l^)-(p^) we 
considered the factor a to have the same value for all the particle species (neutrinos, electrons, positrons). We checked that 
this is a good approximation even for the large asymmetries we are considering. 
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k(e, /x) * 6.5 • 1(T 3 [0.5 (r(&) + rfo,)) + 3(r(-&) + r(-^))] , (21) 
r) ~ 6.5 • 10- 3 [0.5(r(£ M ) + r(&)) + 3(r(-f„) + r(-£ T ))] • (22) 



Numerically, eq. ( ]18[ ) gives: 

MM)-1.45^0"c m ^-^(^) 5 . (23) 

• The expansion of the universe. Oscillations and collisions are ineffective if their scale lengths, l v and l c , are 
larger than the inverse expansion rate of the universe, Ih- In the radiation-dominating regime Ih is expressed 
as: 



1.66,/gT 1 ' 



1.45-10"c„-L(lMfX) 2 , ,24, 

where M p is the Planck mass and g represents the number of relativistic degrees of freedom. We have g = 10.75 
for 1 MeV < T < 100 MeV and g = 3.36 for T < 1 MeV. 



The figure w shows the lengths l v , l c and Ih as functions of the temperature T for ~ 1, r\ e ~ r\ T ~ and various 
values of Am . 




1 o 2 1 o 1 1 1 O 1 o 

T/ MeV 



FIG. 1. The length scales l v , l c and Ih as functions of the temperature T of the electromagnetic radiation in the uni- 
verse. The three thick solid lines represent the vacuum oscillation length l v and correspond, from the upper to the lower, to 
Am 2 = 10" 11 , 10~ r , 10~ 3 eV 2 respectively. The narrow solid line and the dashed line represent the coherence length, l c , for 
the v e — Uft, and the — v T channels respectively. We have taken 7? M ~ 1, rj e — rj T ~ 0. The dotted line represents the inverse 
expansion rate of the universe, Ih- 
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According to the fig. || for Am 2 ~ 10 -7 eV 2 we have l v ~ / c ~ Ih at T ~ 2 MeV. Before this epoch l c & Ih & /«, 
so that collisions are much more efficient than oscillations. The inequality Z c <C implies that the vacuum oscillation 
probability, P c , is suppressed by collisions, as we discussed in this section. As a consequence, for small mixings, 
sin 2 29 -C 1, the flavour composition of the background remains unchanged until T ~ 2 MeV; a partial equilibration 
of the flavours can be realized for large mixings: sin 2 29 <; 0.5. For T < 2 MeV collisions are ineffective, since l c <; Ih, 
and vacuum oscillations develop. 

With the decrease of Am 2 , Am 2 -C 1CP 7 eV 2 , the oscillation length l v increases and, in consequence, for T ^ 2 MeV 
the suppression of oscillations due to collisions is stronger. Even for large mixings the flavour densities are preserved 
until the neutrino decoupling, T ~ 2 MeV. After this epoch oscillations are still suppressed by the expansion rate of 
the universe and become effective, thus changing the flavour composition of the background, only when the oscillation 
length is smaller than the horizon, l v ^ Ih- 

For Am 2 > 10~ 7 eV 2 the inequality l v & l c & Ih is realized before the decoupling. In this circumstance the conver- 
sion probability, P c , is not suppressed, in contrast with the case Am 2 <; 10~ 7 eV 2 . Collisions are effective, thus leading 
to the equilibration of the flavour densities even for small mixing angles. Taking, for instance, Am 2 ~ 10 -3 eV 2 we 
find that equilibration can be achieved for sin 2 29 10~ 3 . Again, after neutrino decoupling the conversion is deter- 
mined by vacuum oscillations. 

For a different choice of the asymmetries at production, e.g. ~ r/ T ~ 1 and rj e ~ 0, the results are similar to those 
in fig. |Tj and we come to analogous conclusions. 

Let us now find the present flavour asymmetries^, 77°, 77^, 77°, for specific neutrino mixings and mass spectra 
motivated by the oscillation interpretation of the solar and atmospheric neutrino anomalies. 
We consider the mixing matrix: 

/ c e -sg \ 
U = s b cq c g ce -s e J , (25) 



.S S Se CgSQ 

where eg = cos 9 and sg = sin 9 and analogous definitions hold for se and cq. The mass eigenstates v\, V2 and v% are 

2% = m) 



related to the flavour ones by the rotation: v a — J2i U a ,i^i- The mass squared differences Am 2 , = m 2 — mf are taken 



t 32 — Am 2 tm ~ 10~ 3 eV 2 and Am^ = Am|, ^ 10 -5 eV 2 according to the currently favoured solutions of 
the solar neutrino problem. Let us first consider Am| x < 10~ 7 eV 2 , as predicted by the LOW and VO solutions. This 
range of Am 21 is the most relevant to the conversion of ultra-energetic neutrinos (see sect. VI). 



According to the results of this section, we identify the following scenario: 

Suppose that before the neutrino decoupling epoch, at T > 2 MeV, a large asymmetry has been produced in one 
flavour while the other asymmetries are initially small: e.g. 77^ = 277 and ?7 e ~ ?7 r ~ 0, with 77 ~ 1. As the universe 
evolves down to T ~ 2 MeV, the muon and tau asymmetries will be equilibrated by the combined effect of oscillations 
and collisions. In the same epochs v e oscillations are still suppressed by collisions. Therefore, the electron neutrino 
asymmetry, 77,= , remains unchanged and, at T ~ Tbbn, we have rj^ rj T = 77 and ?7 e ~ 0. After this epoch neutrinos 
decouple from the thermal bath; collisions become ineffective and the system evolves according to vacuum oscillations. 
During the evolution decoherence occurs due to the spread of the wavepackets. Therefore at the present epoch the 
background neutrinos are in mass eigenstates. With the mixing (|2^) we find the present asymmetries 77^ for these 
states: 

77? ~ ?7sin 2 6 , 77^ ~ 77COS 2 9 , V3-V- (26) 

The corresponding flavour asymmetries equal: 

77 e °^77isin 2 20, 

77° ~ rj (l - ~ sin 2 29 cos 2 G^j , 

77° ~ 77 (l - i sin 2 26* sin 2 6 J . (27) 



4 In the following sections we will consider relatively recent epochs (z ^ 50) for which the flavour composition of the background 
is well approximated by the present one (see sections |ll]| and IV ) . 
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One can see that a large asymmetry is produced in the electron flavour provided that the mixing of the electron 
neutrino is large: sin 2 29 ~ 1. Thus, the electron neutrino asymmetry at the present epoch can be much larger than 
the upper bound (||). Notice also that 77° ~ 77° for 6 ~ 7r/4. 

If equally large asymmetries arc initially produced in the muon and tau flavours, 77^ = r\ T = r\ ~ 1 and ?7 e ~ 0, 
the equality 77^ = r\ T is preserved until the decoupling epoch, T ~ 2 MeV, due to the combined effect of oscillations 
and collisions. The evolution of ?7 e is blocked by collisions for T 2 MeV. After the neutrino decoupling vacuum 
oscillations take place; with the mixing matrix (E5|) we get the same results as in eqs. (p6|)-(p7j). 



1e°^:^?-1- (28) 



For Am 2 ,! — 10 -6 eV 2 and small mixing, sin 2 29 ~ 10 -4 — 10 -3 , according to the SMA solution of the solar neutrino 
problem, the equilibration length l eq is very large. Therefore the electron asymmetry rj e is not equilibrated with the 
muon and tau asymmetries. Again, the present asymmetries are determined by vacuum oscillations which occur after 
the neutrino decoupling and lead to the result (pTj). 

Conversely, for large mixing angle, sin 2 29 ~ 1 and Amfi — 10~ 7 — 1CP 5 eV 2 , as given by part of the LOW solution 
and by the LMA solution regions, equilibration is rapidly realized and one gets: 

2 

3' 

Notice, however, that the results (|2^)-(|2^) depend on the epoch we considered for the production of the large CP- 
asymmetry 77: the equilibration effect of collisions does not take place if the neutrino asymmetries are generated at 
epochs close to the neutrino decoupling epoch, at T ~ 1 MeV. 

2. Evolution in presence of a sterile state 

If a sterile state, v s , is mixed with the three active ones, a general description of the evolution of the neutrino gas 
is complicated and would deserve a detailed study. 

We consider here the specific case in which the sterile neutrino is mixed mainly with one active state only, 
e.g. v ei and the admixture of v s with and v T is negligible. In other words, we consider the mass states 
Vq ~ cos9v e + $m9v s and the orthogonal combination v\ ~ — sin#i/ e +cos#iv Similarly, we take ~ cosc/?i-' AI +sin</jz/ T 
and v-i ~ — sin ipv^ + cos ipv T . This would correspond to v e — v s solution of the solar neutrino problem and — v r 
solution of the atmospheric neutrino anomaly. 

Let us consider the evolution of the v e — v s system. The effective mixing angle in matter, 9 m , can be written as: 

tan2 ^= cos2^ i 2^/A7n^ 

where V is given in eq. ( |l0| ) and E is the average thermal energy of the neutrinos: E ~ a/3T. Numerically, from eqs. 
© and d|) we get: 

sin 29 , . 

tan26» m = -. , 30 

cos 29 - 0.8 • Wa(3Fr, u (T/l MeV) (lO" 3 eV 2 /Am 2 ) 

where we used the expression ti 7 (T) = 2C(3)T 3 /7r 2 for the concentration of photons at the temperature T, with the 
value C(3) — 1.202 for the Riemann zeta function. 

From cq. ( |30| ) it follows that, for r) v 1, T <; Tbbn and Am 2 ^ 1 eV 2 the mixing is strongly suppressed, 
tan2# m <C 1, corresponding to 9 ~ ir/2 (9 ~ 0) if Fr\ v > (Fri^ < 0). Thus, no level crossing is realized before the 
BBN epoch. At T 2 MeV collisions are effective (see fig. [l]) ; however they do not modify ?7 e signifi cantly due to the 
very small value of the mixing and consequently of the conversion probability P c (see section [IB 1, eq. (Uq)). Thus, 



we conclude that no significant flavour conversion occurs and the original value of ?7 e is preserved at least until the 
BBN epoch, even in the case of large vacuum mixing angles. 

As the temperature decreases, T < Tbbn, the mixing angle 9 m approaches rapidly its vacuum value. Taking 77^ ~ 1 , 
F = 2 and Am 2 ~ 10" 3 eV 2 we get tan20 m - tan26» <> 10" 2 for T <, 10 KeV. 



Considering that the propagation of the neutrino states is adiabatic (see section IV B ), we find the present concen- 
trations of v e and v s in terms of the initial density nj^jj 



3 We assume that only active states are initially produced before the BBN epoch, thus n a 
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n? s = n e cos 2 9 , n" e = n e sin 2 9 , if Fr] v > (31) 

n 



- n e sin 2 9 , n° = n e cos 2 (9 , if Fr) v < 0. (32) 



If n e 3> rig, the concentrations of v e can be neglected and relations analogous to (pl|)-(|32|) hold for the CP-asymmetries 
r]° and r\ e . 

The present v^ and v T asymmetries can be found according to the discussion in section [I B 1| . The effect of collisions 
leads to equilibration of Vn and r\ r at T ~ Tbbn'- Vu — Vt = V- At later epochs vacuum oscillations develop, leaving 
this equality unchanged. Thus, we can summarize the present CP-symmetries for the four flavours as follows: 

V° s = rje cos 2 9 , Ve = Ve sin 2 9 , if F<q u > , (33) 

V° s = Ve sin 2 9 , 77° = ?/ e cos 2 9 , if < , (34) 

?7° - ?7° - V ■ (35) 

Having neglected any mixing between v e and the other (active) flavours, we find that the present value of the electron 
neutrino asymmetry is smaller than the one at the BBN epoch, r\ {) e < rj e , thus remaining within the bound given in 



III. HIGH ENERGY NEUTRINO CONVERSION: THE ACTIVE-ACTIVE CASE 

Let us consider three mixed active neutrinos, v e , v^, v r , and find the potential for a beam of high energy neutrinos 
("beam neutrinos") due to the interaction with the relic neutrino background ("background neutrinos"). As discussed 



in the section II B, the flavour composition of the neutrino background changes with time due to the neutrino mixing. 
However, we will focus on neutrinos produced in relatively recent epochs, z 50, when the flavour content of the relic 
neutrino gas has already settled down and does not change with time. 

A. The refraction potential 

According to the section [II B| , decoherence due to the spread of wavepackets implies that the background neutrinos 
are in mass eigenstates. As a consequence, the matrix of potentials, V v , for the beam neutrinos propagating in this 



background is not diagonal in the flavour basis, (v e , v^, v T ). It is possible to check |25| that V v becomes diagonal in 
the basis of the mass eigenstates, [y\, i>2, 1^3), where it can be written as: 

V v = V2G F [(ni - ni) + (n 2 - n- 2 ) + (n 3 - ng)] (36) 

(nif(-s^)-n-J( S ^) \ 

+ V2G F n2/ (_4 2 ))_ n , /( 4 2 )) 

V n 3 /(-4 3) )-ng/(4 J; )/ 

where rii (rq) denotes the concentration of the mass state i>i {vi) in the background and f(s^ ) is the Z-boson 
propagator function: 

f{sf) = . (37) 

(i) 

Here jz and s z are the normalized width of the Z-boson and total energy squared in the Vi — Vi center of mass for 
non-relativistic background neutrinos: 

2Errii _ 2 / E \ / \ T 



sf ~ — ^ ~ 2.4-10 



M 2 7 ' V 10 20 eV 



)(rSr). <38 » 



with rrii being the mass of the neutrino i>{ and E the energy of the beam neutrino. 

The terms in the first line of eq. ([56]) are due to neutral current v — v scattering in the t-channel. The terms in the 
second line of eq. (|36| ) represent the contributions of v. L — z/j scattering with Z-boson exchange in the it-channel, and 
of Vi — Vi annihilation processes. 
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For E 10 20 eV and 1 eV the energy in the i/j — vi center of mass is much below the Z-boson resonance: 



sz 0.03. In this case the propagator function (37) reduces to unity: f(s^) ~ f(~ s ^z) — •"■> ana - ^ nc neutrino- 



neutrino potential ( p6|) becomes energy- independent. 

For extremely high energies, E ~ 10 21 — 10 22 eV, and neutrino mass of order 1 eV the propagator corrections become 
important. However, in this range of energies the absorption effects of the neutrino background are strong |pj,p6|. 
Therefore, the neutrino fluxes at Earth are largely suppressed. In what follows we will concentrate on the low energy 
limit, s% <C 1, which is mainly relevant for applications. 

For a beam of antineutrinos propagating in a neutrino background the potential V p is given by eq. ( |36|) with the 
replacement n% — > nj and vice-versa for all the vi states. 



The fact that the neutrino- neutrino potential matrix, eq. (|36|), is diagonal in the basis of mass eigenstates has a 
straightforward consequence: the effect of refraction consists in a modification of the neutrino effective masses only. In 
terms of the present CP-asymmetries rj? for the mass states Mj of the background we find (for s% <C 1) the following 
corrections: 



IE " 2E " ' " ■ 7V 



Am 



2 



-^1^£ 32 ^^1 + V2G F n 7 (r,° 3 ~ %°) ■ 09) 



As we discussed in section [IB, the present composition of the neutrino background is determined by the initial 



flavour asymmetries, T) e , i]^ and T) T , and by the mixing matrix U of the neutrino system. The expressions of £21 and 



£32 in terms of these quantities can be found from the results of section [IB 2. In particular, with the asymmetries 
@ eq- © gives: 

£si= -2E~ +V *> (40) 

Yji = FjirjV^Gpn^ , (41) 

F 21 = cos 26 , ^32 = sin 2 6 . (42) 

Here we denote as 77 the maximal flavour asymmetry, rj = A/axj^, rj T , r/ e }, which is realized in the background at the 
epoch of nucleosynthesis, T ~ Tbbn', thus 77 is constrained by the bounds (Q). 

In eqs. (^b])-(pT|) the information on the specific mixing matrix and initial composition of the background are 
encoded in the F factors. The dependence of Fj% on the mixing angle in eq. is a consequence of expressing 

the potential (^) is terms of the flavour asymmetry rj, while the background neutrinos are in mass eigenstates. For 
simplicity in what follows we will drop the indexes j,i from the quantities £, V and F in the expressions (|4C|)-([l2|). 



B. The conversion probability 

From the fact that the potential V modifies the effective mass eigenvalues, eqs. (|4C|)-(|l2[), it follows that the 
interaction with the neutrino background does not change the mixing matrix of the neutrino system, which remains 
the same as in vacuum. Conversely, the phase of oscillations is affected by the medium, so that the dynamics of the 
neutrino propagation consists in oscillations with constant depth, given by the vacuum mixing angle, and varying 
oscillation length. The probability P of conversion between two active neutrinos, v a , vp, with mixing angle 9 equals: 

P(t,U) = sin 2 2d sin 2 (^j , (43) 

and the oscillation phase $ is given by: 

= J £(r)dT . (44) 

We denote as U, t the initial and final time of the evolution of the system; £ is given in eq. (f4(i|). Using the scaling 
relations: 

E = E (t Q /t) 2 ^ 3 = E (l + z), V = Vo(t Q /t) 2 = Vb(l + zf , (45) 
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where Eq and Vq are the energy and the potential at the present epoch, z — 0, we get: 



$ = ^vac + $matt , (46) 

with the following expressions for the vacuum oscillation phase, <& vac , and the matter contribution & ma tt'- 



$^0^,^) = — — ^(xs-x?] , (47) 
^mat^a;,^) = Vfao ( — - - J • (48) 



We defined x = t/to and x\ = U/to. 

The matter induced phase, ^ rna tti depends only on the characteristics of the background and on the initial and 
final moments of time. In particular for early production epochs, ^ < 1, one gets: 

$ matt ~ V t — , (49) 

which shows that the phase <5> m att is accumulated mainly at the production time. 

Being independent of Eg/ Am 2 , the phase ^ ma u becomes comparable or even larger than the vacuum oscillation 
phase, $„ ac , at very high energies, Eq/ Am 2 <^ 10 32 eV _1 . Taking Xi = 0.125, corresponding to production at redshift 
z ~ 3, and x = 1 we have $ mc ~ 0.237T for Eo/Am 2 ~ 10 32 eV _1 . This is comparable to the matter phase, 
$ matt ~ 0.297T given by eqs. ([i"8| ) and ( |4l| ) with Ft) ~ 10. As E Q /Am 2 increases the vacuum phase <& vac decreases 
and the total oscillation phase is dominated by the matter contribution § m att- From eqs. (Q) and pq ) we find the 
asymptotic value of the conversion probability^]: 

P{E /Am 2 -> oo) = sin 2 26 sin 2 Q V t (J^ ~ ■ ( 50 ) 

Notice that the expression ( |50|) is insensitive to the change of sign of the potential Vq (i.e. of the product Frf): this 
implies that in the limit of very high energies a beam of neutrinos and one of antincutrinos will experience the same 
matter effect. 

In fig. ^ we show the survival probability, 1 — P, as a function of Eo/Am 2 for neutrinos produced at z = 3 and 
arriving at Earth at the present epoch, with sin 2 29 = 0.5 and various values of the product Frj. The figure was 
produced by averaging the conversion probability, eq. ([l3|), over the interval AEq ~ Eq, keeping in mind the finite 
accuracy in the reconstruction of the neutrino energy in the detector: 



1 



AE 



3Eo/2 



P (Eo) = -r^r / dE'P(E') . (51) 



S /2 



Even though eq. (|5(]) gives a a non-zero value for the conversion probability in the limit Am 2 — » 0, the matter effect we 
are d escribing requires massive non-degenerate neutrinos: Am 2 7^ 0. This condition is necessary for our starting point (section 



[II A, see eq. (Bq)) that the neutrinos in the background are in mass eigenstates (different from the flavour ones), produced 



from flavour states by the spread of the wavepackets during the evolution of the universe. Thus, the expression ([30|) should be 
intended as the high-energy limit of the conversion probability for a given (non-zero) value of Am 2 . 



11 



0.4 I ' ' ' ' 
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Eq/A m 2 (eV" 1 ) 

FIG. 2. The survival probability 1 — P(u a — > vp) as a function of the ratio Eo/Am 2 for various values of Ft). We have taken 
sin 2 28 = 0.5 and production epoch z — 3. 

Let us comment on the figure |[ In absence of asymmetry, Fij = 0, the conversion is given by vacuum oscillations 
(Qmatt — 0). For Eo/Am 2 ^ 10 32 eV -1 the vacuum oscillation phase & vac is very small, thus the conversion 
probability approaches the unity. 

For strongly CP-asymmetric neutrino background the matter induced phase $ r natt is sizeable. In consequence, the 
deviation of the survival probability 1 — P from the value given by vacuum oscillations can be as large as ~ 30%. 
Clearly, as it follows from eq. (|4^) , a strong effect requires a large mixing angle: if sin 2 26 <C 1 the effect of matter on 
the oscillation phase will be unobservable due to the very small amplitude of oscillations. At extremely high energies, 
Eo/Am 2 <; 10 33 eV -1 , the deviation is constant and independent of the sign of Fm according to eq. (|50|). 



IV. HIGH ENERGY NEUTRINO CONVERSION: THE ACTIVE-STERILE CASE 



Let us consider the case in which a sterile neutrino v s is mixed with the active flavours, v a . We discuss here 
two-neutrino mixing; a generalization to a four-neutrino framework will be given in section N% 



A. The refraction potential 



In contrast to the active-active conversion studied in section III , for an active-sterile neutrino system the matrix of 
the refraction potentials is diagonal in the flavour basis (v a , Vs): 



V v = 



v a 





(52) 



In the low energy limit, sz <C 1, the potential V a depends on the flavour asymmetries rj®, if and 77° as follows: 



The potential (|5^ ) can be written in the same general form as eq. (fill): 

V = F<nV2G F n*y , 



(53) 



(54) 
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with the same definition of 77 as the maximal flavour asymmetry in the background at the epoch of nucleosynthesis, 
77 = Max{rj^, rj T , rj e }. The factor F depends on the specific conversion channel and the flavour content of the 
background: 



F = 



(3— e,fi,r 



for v n — v~ channel 



(55) 



Let us consider for instance the conversion of v e to v s . If the v e — v s conversion in the background occurred in the 
resonant channel, the present electron neutrino asymmetry is given by eq. ( |33"| ) . Using this result and eq. ([35]) for 77? 
and 77° we get: 



F = 2 ( 1 + — sin 2 f 
77 



(56) 



where we considered 77^, i] T > rj e . If the v e — v s conversion in the background proceeded in the non-resonant channel Tf e 
is given in eq. ( |34| ) . With this expression one finds a form for the factor F analogous to eq. (^) with the replacement 
sin 2 9 — ► cos 2 9. 

For conversion of antineutrinos the v 
with the replacement F — ► —F. 



v potential has opposite sign: V p = —V v . Thus the expression (p3) holds 



B. The dynamics of neutrino conversion 



Due to the expansion of the universe the cosmological neutrinos experience a potential which changes with time. 
In contrast with the active-active case, the effect of medium changes both the oscillation length and the mixing, eq. 
(^). The dynamics of the flavour transition is determined by the resonance and adiabaticity conditions. 

Consider the resonance condition: 



2EV 1 



= 1 . 



Am 2 cos 29 

Using eq. ( |54j ) and the scaling relations (p5[), from eq. ( |57| ) we get the following relations: 
• The present energy of neutrinos which cross the resonance at the epoch z equals: 

Am 2 \ 10 4 cos26> 



(57) 



En 



10 20 eV 



10 -io eV 2 J Fr)(l + z) 4 



(58) 



For a given Eq and Am 2 the redshift zr at which the resonance condition was realized is given by: 



1 + zr = 10 



cos26> 



F V (10- 10 eV 2 /Am 2 ) (^ /10 20 eV) 



(59) 



• Neutrinos produced at a distance z undergo resonance if their present energy is in the interval: 



E = 10 20 eV 



Am 2 



10~ 10 eV 2 



10 4 cos 26> 
Eq 



1 



{i + zY 



(60) 



Taking Fri = 10, cos26> ~ 1, Am 2 = l(T 10 eV 2 and z = 3 from eq. @ we find E ~ 4 • 10 20 - 10 23 eV. With 
the same values of the parameters and E a = 10 20 eV we get that the resonance condition ( |5q) is satisfied at 
z R ~ 4.6. 



The adiabaticity condition involves the time variation of both the neutrino energy and the concentration of the 
neutrino background. It can be expressed it terms of the adiabaticity parameter at resonance, \Ri as: 
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(Am 2 ) 2 



XR 



AE 



sin 2 26» 





-1 







(61) 



where the subscript "res" indicates that the various quantities are evaluated at resonance, i.e. when the condition 
@ is fulfilled. With the potential Q, using the scalings ( |45| ) and the resonance condition (p7|), we find: 



Xi?, ~ lO -2 ^ tan 2 20(1 + z R )? 



(62) 



For F77 & 10, tan 2 29 = 1 and z# ^ 5 one finds xr & 1-4. Thus, for neutrinos produced at epochs z < 5, we expect 
breaking of the adiabaticity. Notice that \R does not depend explicitly on the neutrino energy and mass squared 
difference; it increases with rj and zr. 

From eq. (p2T) we get the redshift z a corresponding to XR — 2tt ^> 1: 



l + z a 



2?r ■ 10 2 



Frj tan 2 29 



(63) 



If the resonance condition is fulfilled at z > z a the level crossing (resonance) proceeds adiabatically. Taking Fr\ = 10 
and tan 29 = 1 we find z a ~ 15. 

For 77 > 1 and tan 29 < 1 we have z a > z^. Thus, we can define three epochs of neutrino production, corresponding 
to different characters of the evolution of the neutrino beam: 



(i) Earlier epoch: z > z ai when both adiabaticity and the minimum width conditions are satisfied. If also the res- 
onance condition is fulfilled at zr > z a the neutrinos will undergo strong resonance conversion. Otherwise, if the 
resonance condition is not realized (e.g., due to a large value of Am 2 /E) the matter effect can be small. 

(ii) Intermediate epoch: z a > z > Zd- The adiabaticity at resonance is not satisfied (if z a > zr > Zd). At the same 
time the matter width can be large enough to induce significant matter effect. 

Two remarks are in order. (1) The propagation still can be adiabatic in the part of the interval [zd, z a ] outside the 
resonance, and in the whole of it if the resonance condition is never satisfied in this time interval. (2) In monotonously 
varying density the condition for strong matter effect reduces to the adiabaticity condition 0. Therefore, in spite of 
the fulfillment of the minimum width condition, the matter effect can be small for neutrinos produced in the most 
part of the interval [zd, z a \. 

(iii) Later epoch: z < Zd- For neutrinos produced in this epoch the matter effects are expected to be small. 

The fig. H shows the minimum width, resonance and adiabaticity conditions in the z - Frj plane. The minimum 
width condition (|ll| ) is fulfilled in the shadowed region. The lower border of this area corresponds to the curve 
z = Zd(Frj) (eq. Q12|)). For values of Frj and of z in this region one may expect significant matter effect. 

The dashed lines show the values of z and Fr\ for which the resonance condition ( |57| ) is satisfied for neutrinos with 
a given Eo/(Am 2 cos 29) (iso-contours of resonance). 

The solid lines are iso-contours of the adiabaticity parameter: they are contours of constant ratio xrI t arj2 20 
(see eq. (62)). The upper curve corresponds to Xfl/tan 2 20 = 27r, that is, to z = z a for tan 2 29 = 1. For values of 
neutrino production epoch z and Fr] above this contour one would expect resonant adiabatic conversion as dominating 
mechanism of neutrino transformation. For a given Fr] and z the adiabaticity iso-contour gives the value of xrI tan 2 29 
for neutrinos produced at the epoch zi > z and having the resonance at z. In turn, the resonance at z and Frj can 
be satisfied for certain values of Eo/(Am 2 cos 29). It is clear from the figure that strong adiabatic conversion occurs 
for large production epochs, z ^ 10, large asymmetry, Fr\ ^ 10 and large mixing tan 2 29 ~ 0(1)- For Fr\ ~ 2 the 
minimal width condition is fulfilled for large production epochs, z <^ 8, and some effects of adiabatic conversion may 
be seen at z > 15. 
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FIG. 3. The minimum width, resonance and adiabaticity conditions in the z-Frj plane for v a — v a conversion. The solid lines 
are iso-contours of adiabaticity, i.e. of the quantity Xfl/tan 2 26 (numbers on the curves). The dashed lines are iso-contours of 
resonance, i.e. of the ratio Eo/(Am 2 cos 20); the values are given on the curves in units of 10 30 eV _1 . The minimum width 
condition is satisfied in the shadowed region. 

From the above considerations it appears that for realistic parameters a flavour transition of neutrinos occurs either 
due to vacuum oscillations modified by matter effect or by an interplay of oscillations and non-adiabatic conversion. 



C. The conversion probability 

Let us consider neutrinos produced at a given epoch z with a certain flavour v a and propagating in the expanding 
universe with a given constant asymmetry rj. 

We find the v a — v s conversion probability by numerical solution of the evolution equation for two neutrino species 
with the Hamiltonian in the flavour basis (v a , u s ): 

Am 2 nrt0 ()/] , \/-f-\ Am 2 



■^coh26 + V(z) sin 29 

■ sin 26 



H = \ ^ Z L 2 • ^\ ) > (64) 



4E(z) 



where E(z) and V(z) s cale a ccording to eq. (ff5[). 

As discussed in sect. IV B| (see fig. ^), the dynamics of flavour transformation depends on the production epoch 



the resonance epoch, zr, which depends on Eg/Am 2 , and on the value of the adiabaticity parameter at resonance, 
Xr- The figure ^ illustrates the real time evolution of the neutrino states for sin 2 29 — 0.5 and different z, zr, xr, 
which represent different regimes of conversion. 

The solid curve corresponds to production much before the resonance epoch: z > Zr : — 6.8 and weak adiabaticity 
breaking in the resonance, \R — 4-3. The dominating process is the adiabatic conversion which occurs in the resonance 
epoch, tji/to ~ 0.05. The averaged transition probability is close to what one would expect for the pure adiabatic 
case: P a d = 1 — sin 2 9 = 0.85. Weak adiabaticity violation leads to the appearance of oscillations at t > £r. 

The dashed curve corresponds to production close to resonance z <; zr : — 4.2 and strong adiabaticity violation in 
the resonance: \R — 1-2- The dominating process is oscillations in matter with resonance density. At production 
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the mixing is almost maximally enhanced, sin 29 m ~ 1 . The change of matter density leads to slight increase of the 
average conversion probability with respect to sin 2 28 m /2. The decrease of density is fast: the typical scale of density 
change is smaller than the oscillation length, so that maximal depth oscillations do not have time to develop. 

The dotted line shows the same type of regime with stronger adiabaticity violation in resonance. The depth D of 
oscillations is smaller, and the average conversion probability is close to D/2. 
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FIG. 4. The v a — v a conversion probability Pit) as a function of time. We have taken sin 2 29 — 0.5 and three different choices 
of Eo/Am 2 (in units of 10 30 eV _1 ), production epoch z and Fr/. The time t is given in units of the age of the universe, to- 

For further illustration, in fig. ^|we show the evolution in the case of good adiabaticity. Different curves correspond 
to different production epochs: (i) before resonance, z > zr, (ii) at resonance, z = zr, (in) after resonance, z < zr. 
The figures Bfl show similar sets of curves in the cases of moderate and strong violation of adiabaticity. 
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FIG. 5. The v a — v a conversion probability P(t) as a function of time in the regime of good adiabaticity (see fig. |3|). We 
have taken production epochs z earlier, simultaneous and later than the resonance epoch zr. Here sin 2 28 — 0.5, Fr) — 14 and 
Eo/Am 2 = 1.8 ■ 10 28 eV _1 . The time t is given in units of the age of the universe, to- 
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FIG. 7. The same as fig. q for the regime of strong breaking of adiabaticity. Here sin 2 29 — 0.5, Fr/ — 6 and 
E /Am 2 = 4.6 ■ 10 30 eV _1 . 

Let us consider the properties of the conversion probability P(y a — > v s ) for neutrinos produced at epoch z and 
arriving at Earth at the present epoch, z = 0. The probability P depends on z, on the product Fm on the energy and 
mass squared difference in the ratio Eo/ Am 2 , and on the mixing angle 9: P = P(z, Ft], Eq/ Am 2 , 9). As follows from 
figs. the probability is a rapidly oscillating function of z, and also of Eq/ Am 2 . We averaged P over the energy 
resolution interval AEq ~ Eq according to eq. (|5l|). The interpretation of the numerical results can be easily given 
using the z - Frj diagram of fig. ^| 

In fig. ^ we show the dependence of the conversion probability on the production epoch z for different values of Fr] 
and fixed Eq/ Am 2 and sin 2 29. The curves with Frj > represent the resonance channel. For z 1 both vacuum 
oscillations and matter conversion probabilities have oscillating behaviour. For z ^ 2 oscillations are averaged out, so 
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that the vacuum oscillation probability converges to sin 2 29/2 Q A substantial (~ 10%) deviation from the vacuum 
oscillation probability due to matter effect starts at z ~ 1 for Frj ~ 10 and at z ~ 3 for Frj ~ 2. 
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FIG. 8. The v a — v s conversion probability P as a function of the production epoch z for various values of Fr\. From the 
upper to the lower curve: Frj — 20, 10, 6, 2, 0, —2, —6, —10, —20; the dotted line represents the vacuum oscillations probability 
(Fr t = 0). We have taken sin 2 26 = 0.5 and E /Am 2 = 10 31 eV _1 . 

For Frj ~ 6 — 10 and z ~ 4 — 5 neutrinos are produced at densities much higher than the resonance density and 
they cross the resonance at z = 2 — 2.5. The adiabaticity is broken in the resonance, however above the resonance 
the propagation can be adiabatic. For higher asymmetry, Frj <; 10, the adiabaticity starts to be broken near the 
resonance, so that the original flavour state v a ~ v^m will evolve to — {v a + v s )/v / 2. Thus, we have P ~ 1/2. 
With the decrease of z the initial state will deviate from and the conversion probability becomes smaller. With 
the decrease of Frj the adiabaticity starts to be violated earlier (before resonance), so that the transition probability 
decreases. 

For negative values of Frj (or for antineutrinos) the matter effect suppresses the mixing and, in consequence, the 
conversion effect. However, the suppression effect is weaker than the enhancement in the resonant channel. 

Notice that for Frj ~ 10 and z ~ 5 the matter effect can change the vacuum oscillation probability P v by a factor 



For z ~ 2 and Frj ~ 10 the deviation can reach ~ 40% and it equals ~ 20% for Frj ~ 2. 

In fig. H we show the dependence of the survival probability, 1 — P, on Eo/Am 2 for production epoch z = 3, 
sin 2 28 — 0.5 and various values of Frj. Oscillations are averaged for Eo/Am 2 ^ 3 • 10 30 eV _1 ; the averaging 
disappears at Eo/Am 2 ~ 10 32 eV _1 , when the oscillation length approaches the size of the horizon (see also sect. 



Notice that partial averaging exists already at small z due to our integration over AEo- For this reason P does not reach 
its maximal possible value P ma x = sin 2 28. 



of 2: 



{P-P v )/P v ~l . 



(65) 



inland fig. §). 
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FIG. 9. The survival probability 1 — P(v a — v s ) as a function of the ratio Eo/Am 2 for various values of Ft). From the upper 
to the lower curve: Ft) = —20, —10, —6, —2, 0, 2, 6, 10, 20; the dotted line represents the effect of vacuum oscillations (Ft) — 0). 
We have taken sin 2 29 = 0.5 and production epoch z — 3. 
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FIG. 10. The survival probability 1 — P(v a — v s ) as a function of the ratio Eo/Am for various values of the production 
epoch z. We have taken sin 2 29 = 0.5 and Ft) = 6. 



The matter effect increases with Eo/Am 2 . For Eo/Am 2 ^ 5 • 10 30 eV -1 the resonance epoch zr (see eq. (f59|)) is 
earlier than the production epoch of the neutrinos. Thus the neutrinos do not cross the resonance and the matter 
effects is realized mainly in the epoch of neutrino production, when the potential (54) was larger (see fig. |^). The 
value of the effect is determined by the mixing in matter at the production time. With the increase of Eo/Am 2 
the resonance epoch zr approaches the production epoch (see fig. |^). As a consequence the mixing at production, 
and therefore the matter effect, increase. The maximal matter effect is achieved at energies for which the resonance 
condition is fulfilled at production epoch or slightly later (notice that the adiabaticity is strongly broken at resonance) . 
For z ~ 3 this occurs in the interval Eo/Am 2 ~ 10 31 — 10 32 eV -1 . For z ~ 5 maximal matter effect is realized at 
Eo/Am 2 ~ (5 - 7) • 10 30 eV" 1 (fig. 0). 
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In fig. |ll| we show the dependence of the matter effect, i.e. the difference P — P v , on the quantity Frj for various 
values of the mixing angle. For the parameters used in the plot the neutrinos are produced close to the resonance 
and the adiabaticity is strongly violated in the resonance. The matter effect can be estimated as the deviation of the 
jump probability from 1: 

1-P LZ ~1- cxp(-^ X fl/2) • (66) 
In our case \R <C 1, so that the matter effect is proportional to Fr]: 

P-P V ~^ X R^ Ft} tan 2 29 , (67) 

according to eq. (|62|). This explains the linear increase of the matter effect with r\ and F. 
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FIG. 11. The deviation with respect to the vacuum oscillation probability, P(y a — v a ) — P v as a function of the product Fr] 
for various values of sin 2 28. We have taken production epoch z = 3 and Eq /Am 2 — 10 31 eV _1 . 

In fig. [l^ we show the dependence of P — P v on the mixing parameter sin 2 29 for different values of the ratio 
£o/Am 2 and fixed production epoch z = 3 and Fr] = 6. The neutrinos are produced in the resonance epoch or after 
it depending on their energy. For small mixing the matter effect is proportional to the mixing parameter sin 2 29 m at 
the production time. This explains the linear increase of the effect with sin 29 (sin 2 29 m oc sin 29) and with Eq / Am 2 
(for Eo/Am 2 ~ 10 31 eV _1 the production epoch coincides with the resonance one). For maximal mixing, sin 2 29 = 1, 
the average probability takes the value P = 1/2 independently on adiabaticity violation [p7| . Therefore in this case 
P — P v = 0. The maximum deviation from vacuum oscillation effect is realized at sin 2 29 ~ 0.65. 
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ne 12. The deviation with respect to the vacuum oscillation probability, P{y a — u s ) — P v as a function of sin 2 29, for various 
values of Eg/Am 2 (in units of 10 30 eV _1 ). We have taken production epoch z = 3 and Frj = 6. 



V. CONVERSION EFFECTS ON DIFFUSE NEUTRINO FLUXES 



The results we have discussed in the sections III B and [V C describe the conversion effect for a beam of neutrinos 



produced by a single source at a certain epoch z. Presently, the possibilities of detection of neutrinos from single 
sources are limited to objects with redshift z <C 1. For these neutrinos no substantial effect is expected^. There is a 
hope, however, to detect the diffuse (integrated) neutrino flux which is produced by all the cosmological sources. For 
this flux matter effects can be observable. 

In what follows we will calculate the ratio F a (Eq) / Fq (E ) , where F a (E ) and Fq(E ) are the present diffuse fluxes 
of neutrinos of given flavour, v a , and a given energy, Eq, with and without conversion. The ratio can be written as: 



F a (E ) 
F a a (Eo) 



l-P a {E ) , 



where P a is the averaged transition probability: 

P Q (Eo) ' 



Fq(E ) Jo 



dF£(E ,z) 
dz 



P a (E ,z)dz . 



(68) 



(69) 



Here P a ( Eo, z) is the transition probability for neutrinos produced in the epoch z, which has been discussed in sections 
[II B and IV C| . The quantity (1Fq(Eq, z) is the contribution of the neutrinos v a produced in the interval [z, z + dz] 
to the present flux in absence of oscillations. 

We first derive the general expression for the differential flux dF^(_Eo, z). Let f(E) be the flux of neutrinos generated 
by a single source. Then the total number of neutrinos produced in the unit volume in the time interval [t, t + dt] 
with energy in the interval [E, E + dE] can be written as: 



f(E)n(t)dEdt , 



(70) 



where n(t) is the concentration of sources in the epoch t. The contribution of these neutrinos to the present flux 
equals: 



dF a (E Q , z) = ±f(E)n(t)(l + z)- 3 ^-dt 



(71) 



Some effect can appear due to conversion in halos of galaxies and of clusters of galaxies 
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where c is the speed of light and the factor (1 + z) 3 accounts for the expanding volume of the universe. Transferring 
from t to z- variable we get: 

dF«(E , z) = ^f(E)n(z)(l + z)-^ 2 ^-dz . (72) 

The relation between the energy E and the present neutrino energy E includes, in general, effects of energy losses 
and of redshift. Neglecting absorption we have dE/dEo — (1 + z). The density of sources, n(z), can be expressed in 
terms of the comoving density n c as n(z) = (1 + z) 3 n c (z). Notice that n c — const, if the number of sources in the 
universe is constant in time. Thus, the evolution of sources is described by the dependence of n c on the redshift z. 
In terms of n c and Eq we get finally: 

dF£(E , z) = ^/(£ (1 + z))n c {z){\ + z)-^ 2 dz . (73) 
Inserting dF a (Eo, z) in eq. j69| ) we find: 

Pa{Eo) =FfW)^ I f^ 1 + ^))n c {z){l + z)-^P a {E Q ,z)dz , (74) 

and Fq(Eq) is given by the same expression with P a = 1. 

In what follows we will calculate the survival probability 1 — P a for various possible sources of high-energy neutrinos, 
assuming certain forms for the produced flux f(E) and the concentration of sources n c . 



A. Conversion of neutrinos from AGN and GRBs 



There is an evidence that cosmological sources like GRBs and AGN were more numerous in the past. In particular, 
the density of GRBs evolved as [M: 



r (i + z) 3 z < Zp 

n c (z) oc < (1 + z p ) 3 z p < z < z max 

V ~ Z > Zmax 



(75) 



where z p is estimated to be z p ~ 1 — 2 [^9| . The energy spectrum of neutrinos from GRBs scales as a power law |30| 



/(£)« 



E 2 E 2 {l + zf 



(76) 



Combining eqs. (|7q) and (|76[) with (74) we find the averaged probability: 



(l + z)- 1 / 2 P a {E 0l z)dz+ / (1 + z)-^ 2 P a (E ,z)dz 



(77) 



where the normalization factor N p is given by the expression in square brackets with P a — 1. According to eq. ( |77| ) 
the contribution of the recent epochs to the present flux is enhanced in spite of the the larger number of sources in 
the past. This leads to suppression of the matter effects, which are more important at large z. 

The figure [l3] shows the averaged survival probability , 1 — P a , for v a — v s conversion channel, as a function of 
Eo/Am 2 for different values of Frj. We have taken z p = 2 and z max = 5. The averaged probability is rather close to 
the non-averaged one (see fig. ||) for neutrinos produced at z ~ Zp = 2. Indeed, the contribution to the flux from the 
earlier epochs, z <; z p is strongly suppressed, according to eq. ([770 ■ The integration over z leads to some smoothing 
of the oscillatory behaviour of the probability. The deviation of the ratio F a (Eo)/ Fq(Eq) from its vacuum oscillation 
value can reach ~ 25%. Maximal effect is realized for Frj ~ 20 in the resonance interval Eq / Am 2 ~ (1 — 5) • 10 31 eV _1 . 
For F t] ~ 2 the effect is about (3 - 4)%. 
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FIG. 13. The averaged survival probability for v a — v s channel, 1 — P a , as a function of the ratio Eo/Am for the dif- 
fuse flux of neutrinos from GRBs. The curves correspond to various values of Frj. From the upper to the lower curve: 
Frj — — 20, —10, —6, —2, 0, 2, 6, 10, 20; the dotted line represents the effect of vacuum oscillations (Frj = 0). We have taken 
sin 2 26» = 0.5. 

For conversion between active flavours the results are shown in fig. the deviation of the survival probability 
from the value given by vacuum oscillations can be as large as ~ 10% for large asymmetry, Frj <^ 10, and high energies, 
Eo/Am 2 ^ 10 32 eV -1 , for which the matter- induced oscillation phase $ ma tt dominates over the vacuum oscillation 
phase <&„ ac . 
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FIG. 14. The averaged survival probability for v a — vp oscillations, 1 — P a , as a function of the ratio Eo/Am 2 for the diffuse 
flux of neutrinos from GRBs. The curves correspond to various values of Ft]. We have taken sin 2 29 = 0.5. 

The astrophysical data about AGN indicate that the distribution of these objects has maximum at z ~ 2 fl3"l| , 
with a rapid decrease of the concentration with z. The power law f(E) cx E~ 2 is a good approximation for the most 
energetic part of the spectrum p2| . For these reasons, in the case of AGN the results are similar to those discussed 
here for neutrinos from GRBs. 
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B. Conversion of neutrinos from heavy particle decay 



Very heavy particles, with mass of the order of the grand unification scale, are supposed to be produced in the 
universe by topological defects, e.g. in monopole-antimonopole annihilation, cosmic strings evaporation, etc. 
These particles would then decay very quickly, with lifetime r <C to, into leptons and hadrons. Neutrinos may be 
produced cither directly, as primary decay products, and/or as secondary products from decays of hadrons. 

Let us calculate the contribution of the neutrinos produced in the epoch z to the present flux: dF^{Eo, z)/dz. In 
assumption of very fast decay of the heavy particle, X (so that the production epochs of X and of the neutrinos 
coincide), we can write the total number of neutrinos produced in the unit volume in the time interval [t, t + dt] with 
energy in the interval [E, E + dE] as: 

dn x (t) dN 

-^lLE dtdE ' (78) 

where dnxit) is the number of X particles produced in the interval [t, t + dt] in the unit volume, and dN v is the 
number of neutrinos in the energy interval [E,E + dE] produced by a single particle X. The contribution of the 
neutrinos produced in the epoch t, eq. (|78|), to the present neutrino flux is: 

, c dnxit) dN v , , „ dE , 

dF ^^ = ^^r^ 1+z ^ dE dt > < 79 > 

where we have taken into account the expansion of the universe. In terms of the redshift z we get: 

dF "(E ,z) = ^^( Z )^(E (1 + z))(l + z)-^dz , (80) 

where we used also the relation E = Eo(l + z). 

The production rate of the X particles can be written as: 

dn x (t) 4+ 6 _- 



dt 



oc t- i+p oc (1 + zf-iP , (81) 



where p = 1 for monopole-antimonopole annihilation and cosmic strings |33|] and p = 2 for constant comoving 
production rate. 

For the fragmentation function of neutrinos we take a power law: 

^<xE«=ES(l + zr. (82) 
If the neutrinos are produced mainly by hadronic decays the fragmentation function has a polynomial form p4| . The 



leading term of the polynome gives the expression (82) with a = —3/2. 



Inserting the expressions from (BQ), (|8l|) and (B2h in eq. (p9) we get: 



and, for p — 1 and a = —3/2: 



P a (E ) = ^- /(I + :!"--"-"/>„(£;, -.),/: . (,s:j) 



P a {E Q ) = ^- III + :) 7 />.(£■„.:),/: . (M) 



Here N p is a normalization factor. 



We perform the integration ( |84| ) starting from the absorption epoch z a b s . The contribution of the neutrino flux 
produced at z z a b s is very small due to absorption^. The dominant absorption processes are v — v and v — v 
interaction with the neutrino background. The absorption epoch is given by: 



9 Clearly, the energy of the neutrinos at production can not exceed the mass of the parent particle, X. This gives a 
further constraint on the upper integration limit: 1 + z max ^ mx/Eo- Stronger bounds can be found in some specific 
production mechanisms: taking, for instance, X — > tv + tv~ and subsequent production of neutrinos by pion decay, one gets 
1 + Zmax ^ 0.2134mx/£o |35| ■ For Eo ^ 10 22 eV and mx ~ 10 16 GeV this gives the constraint z ma x 200, which is weaker 
than the one given by absorption. 
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1 + Zabs = 



dabs 

du 



(85) 



where djj is given in eq. (Q) and d a b s is the the absorption width, which depends on the v — v energy squared in the 
center of mass, sz, (see eq. (^8h). Taking, for instance, E ^ 10 22 eV and m v 0.05 eV, we have sz 0.1, and the 
corresponding absorption width is d a b s <; 1.5 • 10 34 cm -2 p^j. With this value and 77 ~ 10 eq. ( |85| ) gives z a b s ~ 50. 

In the figure ^| we show the averaged survival probability 1 — P a for v a — v s conversion channel, as a function of 
Eq/ Am 2 for different values of Fij. One can see that, in contrast with the case of neutrinos from GRBs, the deviation 
of the ratio F a (Eq) / Fq (Eq) from the value given by vacuum oscillation is significant (larger than ~ 10%) in a wide 
range of energies: E /Am 2 ~ 10 26 — 10 32 eV -1 . 




FIG. 15. The averaged survival probability for v a — v s channel, 1 — P a , as a function of the ratio Eo/Am 2 for neutri- 
nos from the decay of heavy relics. The curves correspond to various values of Frj. From the upper to the lower curve: 
Frj — —20,-10,-6,-2,0,2,6,10,20; the dotted line represents the effect of vacuum oscillations (Frj = 0). We have taken 
sin 2 26» = 0.5. 



For a given value of Eo/Am 2 the matter effects are determined by the corresponding resonance epoch, zr, and 
adiabaticity in resonance. For Eq/ Am 2 < 10 29 eV -1 the resonance was realized at zr > 10, when the adiabaticity 
condition was fulfilled (see fig. ||). Therefore, the matter effects are do mina ted by resonant adiabatic conversion 
which occurs for neutrinos produced at z > zr ~ 10. As discussed in sect. [VC, these neutrinos undergo almost total 
conversion (see fig. |]), however, their contribution to P a is suppressed according to eq. (|84|). 

For Eq/ Am 2 *J 10 29 eV -1 the resonance is realized at zr ^ 10, when the adiabaticity is broken (fig. ||), so that 
the matter effect is mostly due to non-adiabatic conversion and oscillations in the production epoch. 

The maximal effect is realized in the interval Eq/ 'Am 2 ~ 10 29 — 5 • 10 31 eV -1 ; the relative deviation of 
F a (Eq) / Fq (Eq) with respect to the vacuum oscillations value equals ~ 10% for Fr\ = 2 and can be as large as 
50% for Frj = 20. 

The figure |l6| shows the average survival probability, 1 — P a , for active-active conversion. We see that, similarly to 
what discussed for neutrino from GRBs, a substantial (~ 15%) matter effect requires large asymmetry, Frj <; 10, and 
very high energies, Eq/ Am 2 <; 10 32 eV -1 , for which the matter contribution to the oscillation phase is dominant. 
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FIG. 16. The averaged survival probability for v a — up oscillations, 1 — P a , as a function of the ratio Eo/Am 2 for the diffuse 
flux of neutrinos from the decay of heavy relics. The curves correspond to various values of Frj. We have taken sin 2 26 = 0.5. 

VI. OBSERVABLE EFFECTS 

Let us consider the experimental signatures of matter effects on neutrino propagation. The observable effects depend 
on the specific scheme of neutrino masses and mixings and on the initial flavour composition of the neutrino flux. 



A. Conversion of cosmic neutrinos and neutrino mass schemes 



As follows from the analysis of sections VA-VB, a significant matter effect on active-active oscillations of high- 
energy neutrinos requires: 



E 



Z 10 32 eV -1 



(86) 



This is the condition for which the matter-induced oscillation phase, $ ma tt, dominates over the vacuum one, <& vac 
(see section IIIB). For conversion into a sterile neutrino the matter effect is substantial in the ranges: 



E ( 10 30 eV" 1 for AGN, GRBs 
Am 2 \ 10 28 eV -1 for heavy relics decay 



(87) 



For E ^ 10 21 eV, the conditions (86)-(B7ft imply 



Am 2 £ i 10 ~" eV " for - "f> 



ii eV 2 
10- 7 eV 2 for v a - v s 

For both the active-active and active-sterile channels the mixing angle should be large enough and, for v a — v 8 , not 
too close to maximal (see fig. |l2"|): 



0.1 <, sin 2 26> < 0.95 . 



(89) 



In the three neutrino schemes which explain the solar and atmospheric neutrino data the effect can be realized for 
Ve-Vr/v^ mixing and the vacuum oscillation (VO) solution of the solar neutrino problem. If the LMA, the SMA or 
the LOW solution are confirmed, the effect of medium on vacuum oscillations of cosmic neutrinos can be neglected. 

In presence of a sterile neutrino the conditions ( p7[ ) - ( |89| ) can be realized in a number of situations. Oscillations 
of electron neutrino into a sterile state, v e — v s , with Am? <J 10~ n eV 2 and mixing close to maximal represent a 
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possible solution of the solar neutrino problem p&. Another possibility is to consider, e.g., the hierarchical mass 



spectrum with 7713 ~ \J Am^ fm , mi ~ ^/Atoq, mi ^ 10 4 eV and mo < mi so that Am^ ~ m\ ^ 10 7 cV 2 . In 

the simplest case the sterile state is mixed only in the lightest mass eigenstates v\ and v$. The mixing angle is only 
weakly restricted by the solar neutrino data|^. 

B. Flavour composition of detected fluxes 

Let us consider the numbers of events N a and N® induced in a detector by neutrinos of different flavours a with 
and without conversion respectively. These quantities are determined by the present fluxes F a and Fq (see section 

0) ; if the detector provides total energy reconstruction and optimal event selection, the flavour composition of the 
numbers of detected events coincides with that of the fluxes. 

In what follows we consider two possible types of flavour composition for the numbers of events in absence of 
conversion: 

1) . CP-symmetric: N° = {a = e,fi, t). As far as flavour content is concerned we take the normalized numbers of 
events: 

(N° e ,N°,N?) = (1,2,0) 

(NlN°,N°) = (1,2,0) . (90) 

Such a flavour composition is expected for neutrinos produced by the decays of n + and tt~ mesons, which in turn 
appear in the process X — > ir + ir~ (see section \ 



2). CP-asymmetric: iV° ^ iV?. We consider 

(Ar°,jV°,AT°) = (1,1,0) 

(Ni,N°,N°) = (0,1,0). (91) 

This flavour composition is realized for neutrinos produced by the scattering of highly energetic protons on a photon 
background, where the ir + decay gives the dominant contribution. The p-j interaction is supposed to be the main 
mechanism of neutrino production in GRBs p0| . 

Neutrinos of different flavours produced in the same decay reaction (X or it decay) share the energy of the parent 
particle equally with good approximation. Therefore the produced fluxes of neutrinos and antineutrinos of different 
flavours have the same energy dependence, and, in absence of conversion, the ratios N® /N®, N® /N®, N®/N® are 
expected to be energy-independent. 

In presence of vacuum oscillations the ratios of numbers of events are approximately independent of energy in two 
intervals: 

(i) E /Am 2 <, 5 ■ 10 30 eV" 1 (see figs. where oscillations are averaged out. 



(ii) Eo/Am 2 ^ 5 • 10 32 eV -1 , where the vacuum oscillation phase is very small or, equivalently, the vacuum oscillation 
length exceeds the size of the horizon. In this case the conversion probability is negligibly small and the ratios of 
numbers of events approach their values in absence of oscillations. 

The effects of vac uum osci llations are modified by the interaction with the neutrino background. According to the 



results of sections VA-VB wc find that 



1. The energy dependence of the ratios of numbers of events in the interval (i) would be a signal of active-sterile 
conversion with matter effects^. 



10 No restriction exists for Am 2 <C 10~ n eV 2 ; bounds follow from the solar neutrino data for Amf <; 10 -11 eV 2 (in this 
case the solar neutrino data should be treated in three neutrino context). 

xl If some difference exists in the energy dependences of the original fluxes of neutrinos of different flavours, this would appear 
in the total number of events iV t ° ot = N° , in contrast with the effect of neutrino conversion. Thus, the energy-dependence 
of ratios of numbers of events due to matter effects can be distinguished. 
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2. The deviation of the ratios of numbers of events in the interval (ii) from the values expected in absence of 
conversion would indicate matter-affected active-active oscillations. Two elements, however, will make the 
identification of the effect difficult: its appearance at very high energies, close to the end of the predicted 
spectra of ultra-high energy neutrinos, and the uncertainties on the flavour composition of the neutrino fluxes 
at production. 

We notice an interesting aspect: the interaction with the neutrino background produces strongly different effects 
on active-active and active-sterile oscillations. Thus the observation of such effects would neatly distinguish between 
the two channels. In particular, the observation of the characteristics described in 1. would give indication of the 
existence of a sterile neutrino. 



C. Ratios of numbers of events: active neutrino mixing 

For three neutrino flavours, v e , v T , the relation between the numbers of events N a and iV° can be expressed as: 

Nu = VN° , (92) 

where: 

N? = (N e ,N°,N?) 

N v = (N e , Njj,, N T ) , (93) 
and V is the matrix of conversion probabilities: P a [3 = P(y a — > vg), (a, = e, /a, t) 



As an example we consider the scenario introduced in section [I B 1 , in which the solar neutrino problem is solved by 
Ve — v^lvr vacuum oscillations with Atoq = Am^ ~ 10~ n eV^ and the atmospheric neutrino anomaly is explained 
by v y, ~ v t oscillations with Am 2 tm — Am§ 2 ~ 10~ 3 eV 2 . The mixing matrix is given by eq. (|2^). 

Since the values of Am§ 2 and Am 2 ^ (Am 2 ^ ~ Am 2 2 ) are out of the range of sensitivity to matter effects (see sect. 



VIA) the the oscillations due to Am^ and Am| 2 are described by the average vacuum oscillation probability. The 
neutrino background influences the v\ — v-i system only. In these specific circumstances matter effects show up in 
the conversion of v e = cos2#i/i + sin2#z/ 2 into the orthogonal state v' = — sin 26*1^ + cos26>i/ 2 . We denote by P the 
corresponding two-neutrino conversion probability. Taking the maximal mixing O = 7r/4 in the matrix ( p5| ) we find 
the conversion matrix: 

P/2 P/2 \ 

V= \ P/2 1/2 -P/4 1/2-P/4) , (94) 




1/2- P/ A 1/2 -P/4, 

and an analogous expression for the matrix of probabilities for antineutrinos with the replacement P — > P, where P 
represents the D e — > v' conversion probability. 

Taking the CP-symmetric flavour composition (po|), from eqs. ( |94| ) and ( p2|) we find that the conversion probability 
P cancels in the expression of the numbers of events, N a . Equal numbers of events for the three flavours are predicted 
independently of matter effects: N u = N„ = (1, 1, 1). 

For the CP-asymmetric composition (|9l]) we obtain: 

N v = (1 - P/2, 1/2 + P/4, 1/2 + P/4) 

N„ = (P/2, 1/2 - P/4, 1/2 - P/4) . (95) 

Since the present detectors do not distinguish neutrinos from antineutrinos, we consider the sums of the events induced 
by v and v. From eqs. ( j95[) we find: 

N v + N v = (1 - P/2 + P/2, 1 + P/4 - P/4, 1 + P/4 - P/4) . (96) 

Two comments are in order. First, equal numbers of events induced by the muon and tau neutrinos are expected, 
with no dependence of ratios on matter effects: (N^ + Np)/(N r + Nf) — 1. Conversely, matter effects are present 
in ratios involving the electron neutrino. Second, if P = P the conversion probability cancels in (^6|) and one gets 
N v + N„ = (1, 1, 1). This circumstance is realized in absence of matter effects (Frj = 0) or in the extremely high 
energy limit, Eo/Am 2 ^ 10 33 eV -1 , in which the asymptotic value (J50|) for the conversion probability is realized (see 
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also fig. Therefore, the matter effect could be revealed by a deviation from the equality of number of events for the 
three flavours in the narrow interval Eq/ Am 2 ~ 10 32 — 10 33 eV _1 in which P and P are unequal and have significant 
deviation from the vacuum oscillation probability. 

Considering, for instance, the ratio of e-like over non e-like events we find: 

R = *«±_* = l-P/2 + P/2 

Nft + Nfi + N T + Nf 2 + P/2 - P/2 y ' 

The deviation of R from its value R p = 1/2 without oscillations is entirely due to matter effects and equals: 

- -§A , (98) 



where A = (P — P)/2. The relative deviation (|9§|) amounts to ~ 15% for A ~ 0.1. Similar conclusions are obtained 
for other ratios of numbers of events. 

Results arc different if the mixing angle O in the matrix (p5|) is not maximal. In the extreme case = the 
problem reduces to two-neutrino conversion. In the limit Eo/Am 2 <; 10 33 eV _1 , for both the compositions ( |90| ) and 
© we get: 

R - R p „ j P m 

R p ~ 2 ' 1 > 

where we considered P ~ P. Taking P ~ 0.1 the deviation (|99| ) equals ~ 15%. 



D. Extension to four neutrinos 



An example of four neutrino scheme with sterile neutrino, v. 
present in the two light mass eigenstates, vq and i>i, so that in the bases v, 
mixing matrix takes the form: 



was introduced in section |VI A 
(v s ,v t 



Vr), 



u° = 



Ccj, 

-ses^j\[2 
\-ses$/V2 



S<f, 

seep/ \[2 
sec^/ V2 









so 
c e /V2 -1/V2 
c e /V2 1/V2 J 



the sterile state is 
= (^0,^1,^2,^3) the 



(100) 



= COS0, 

s e Vu,/V2 + 



SllK 



co 



cos 9, sg = sin 9. The angle <f> describes the mixing between v s and the state 
sqv t I 'y/2. Analogously to the previous case, we consider Am\ a ^ 10 -7 eV 2 and all the other 



where 

V = CgV e - 

mass splittings to be much larger than this value, so that the interaction with the neutrino background affects the 
propagation of the vq — v\ system only. As a consequence, the matter effect modifies the angle <f> only; the changes of 
are negligibly small. Again, the dynamics of the four neutrino system is reduced to the evolution of the two states 
v s and v. Introducing the conversion probability P = P(y s — > v), we find the matrix of probabilities (see eq. (92)): 



V 



fl-P 

2p 



s\P/2 
sjP/2 



At 



sic 



F4(l-P) 

1(1 -P/2) 
2 (1 - P/2) 



s 2 e P/2 
14(1 -P/2) 



(1 



P)]/4 
P)]/4 



slP/2 
S 2 c 2 (l-P/2) 
[l + 4 + sj(l-P)}/4 
[l + 4 + sj(l-P)}/4 



(101) 



Taking the CP-symmetric composition ( p0| ) and assuming that no sterile neutrinos are produced, N® = 0, from eq. 
(101) and (pi) one gets the numbers of events: 



N v = (P, 1 - c 2 P, 1 - ,s 2 P/2, 1 - ,s 2 P/2) 
N P = (P, 1 - c 2 P, 1 - s 2 P/2, 1 - s 2 P/2) 



(102) 



As in the three neutrino case, we have (iV M + Np)/(N T + Nf) = 1 independently on matter effects. Notice that in 
the total numbers of events 7V„ + N„ the conversion probabilities appear in the combination P + P: since the matter 
effects have opposite signs for neutrinos and antineutrinos, they partially cancel in this quantity. 

Introducing the deviation from the averaged vacuum oscillation probability, Sp = P + P — 2P V , we compute the 
ratio: 
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R : 



N e +N s = l-4(P v + S P /2) 

N, t + N fi + N T + N t 2-s 2 (P v + 6 P /2) 



The relative deviation of this ratio from the value given by vacuum oscillations equals: 

R-Rv s P 4-4/ 2 



R r 



2 (l~s 2 P v /2)(l-c 2 P v ) 



(103) 



(104) 



Taking 5 P ~ 0.1, sj 
small 9: ci ~ 1, si i 



1/2 and P v ~ 0.4, eq. (104) gives a deviation of ~ 2%; the effect is larger, ~ 10%, for 



For the CP-asymmetric composition (|9l]) we get: 

Af° = (P(c 2 + 4/2), 1 - 44 - c 2 e P(c 2 e + 4/2) , 

(1 + s 2 c 2 )/2 - s 2 P(c 2 + s 2 /2)/2, (1 + s 2 c 2 )/2 - s 2 P(c 2 + s 2 /2)/2) 
N° = ( S 2 P/2,s 2 c 2 (l~P/2), 

[1 + 4 + 4(1 - P)]/4, [1 + 4 + 4(1 - P)]/4) . 

For the ratio, P, of the e-like over non e-like events one gets: 



R — R v 

Py 



[5(1 - 4/2) + 54/2] 



c 2-4/2 



(l- s 2 P„/2)(l-c 2 P„) 



(105) 



(106) 



where 5 = P — P v and 5 = P — P v . 

with the values S ~ 0.1, 5 ~ —0.05, P v ~ 0.4 and small mixing, Cg ~ 1, s| — 0> the deviation ( |106| ) equals 
15% similarly to the case of CP-symmetric composition, eq. (|104|). The effect is smaller, ~ 2%, for large mixing, 



s 2 ^c 2 ^ 1/2. 



Our estimation, 10 — 15% effect, gives some hope that the discussed phenomenon will be observed in future large 
scale experiments with event rates ~ 1000 events/year. 



VII. CONCLUSIONS 



We have studied matter effects on oscillations of high energy cosmic neutrinos. The only known component of the 
intergalactic medium which can contribute to such an effect is the relic neutrino background provided that it has large 
CP (lepton) asymmetry. 

The mixing modifies the flavour composition of the relic neutrino background. Considering atmospheric and solar 
neutrino- motivated mixings and mass squared differences we find that, if large asymmetries in the muon and/or tau 
flavours are produced before the BBN epoch, they are equilibrated by the combined effect of oscillations and inelas- 
tic collisions, so that rj^ ~ r/ T . The asymmetry in the electron flavour, rj e , can be equilibrated with r]^ and r\ T for 
Am 2 ^ 10~ 7 eV 2 . For Am 2 & 10~ 7 eV 2 v e — v^/vt oscillations are suppressed by collisions and/or by the expansion 
rate of the universe, thus leaving r] e unchanged, at least until the BBN epoch. At later epochs oscillations develop 
and large asymmetries in the muon and/or tau flavours can be efficiently converted into u e asymmetry. Therefore at 
present the values of the asymmetries for the three flavours can be comparable. This allows one to reconcile possible 
large lepton asymmetry in the v e flavour at present, r\ e ~ 1, with strong constraint on rj e from nucleosynthesis. Active- 
sterile conversion is ineffective until the BBN or later, due the matter-induced suppression of the v a — v s mixing. 

The dynamics of high-energy neutrino conversion in the CP-asymmetric neutrino background has been considered. 
For conversion between active neutrinos the matter effects consist in a modification of the vacuum oscillation length. 
The effect is significant for large mixing angle, sin 2 20 ^ 0.3, and high energies, Eq/ Am 2 <; 5 • 10 32 eV -1 , for which 
the matter contribution to the oscillation phase dominates over the vacuum oscillation one. In these circumstances 
the conversion probability can differ by ~ 30% from the vacuum oscillations value. 

For active-sterile conversion the matter effects can be important in the interval Eq/ Am 2 ^ 10 28 — 10 32 eV -1 , for which 
the resonance condition is satisfied. For the majority of realistic situations (with z ^ 10, rj ^ 10), the adiabaticity 
condition is broken. This implies that the matter effect is reduced to non-adiabatic level crossing or enhancement 
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(suppression) of mixing and therefore of the depth of oscillations in the production epoch. The relative change of 
the conversion probability can be as large as 20 — 50%. For extreme values of the asymmetry and large redshift of 
production the adiabatic conversion can take place with almost maximal conversion probability. 

We calculated the effect of conversion on the diffuse fluxes of neutrinos produced by GRBs, AGN and the decay of 
super-heavy relics. For neutrinos from GRBs and AGN the relative deviation of the flux due to matter effects with 
respect to vacuum oscillations can reach 20%. For neutrinos from heavy particle decay the effect can be larger: up to 
~ 40%. 

Possible signatures of matter effects consist in the deviation of ratios of numbers of observed events, 
Ne/N^, N e /N T , N^/N r , from the values predicted by pure vacuum oscillations. Presumably, neutrino mixings and 
masses will be measured in laboratory experiments and vacuum oscillations effects will be reliably predicted. 

For conversion into a sterile state one expects also a characteristic energy dependence of the ratios which in principle 
will allow to distinguish matter effects from the uncertainties in the flavour content of original neutrino fluxes. 

For illustration purpose we estimated observable effects for two possible schemes of neutrino masses and two dif- 
ferent flavour compositions of the detected fluxes in absence of conversion. In a scheme with three flavour states only 
and parameters in the region of VO solution of the solar neutrino problem we found that the deviation of ratios of 
numbers of events from their vacuum oscillation values can be of ~ 10%. Similar conclusion is obtained for schemes 
with an additional sterile neutrino. 

Clearly, more work is needed to clarify the possibilities to observe the effect under consideration. In any case, new 
large scale detectors with relatively high statistics (<~ 1000 events/year) are required. 

The detection of matter effects on fluxes of high-energy neutrinos would be an evidence of large CP (lepton)- 
asymmetry in the universe. As follows from our analysis, asymmetries of order r\ ~ 1 can be probed in these studies. 
Clearly, the observation of such a large asymmetry will have far going consequences for our understanding of the 
evolution of the universe. 



[1] see e.g. the reviews: R.J. Protheroe, Nucl. Phys. Proc. Suppl. 77 (1999) 465. [astro-ph/9809144], E. Waxman, hep- 

ph/0009152, and references therein. 
[2] For a review, see P. Bhattacharjee, in College Park 1997, Observing giant cosmic ray air showers, 168-195; Edited by J.F. 

Krizmanic, J.R. Ormes, R.E. Streitmatter. Woodbury, AIP, 1998. 536p. (1997), [astro-ph/9803029]. 
[3] T.J. Weiler, Astropart. Phys. 11 (1999) 303, [hep-ph/9710431]; D. Fargion, B. Mele and A. Salis, Astrophys. J. 517 (1999) 

725, [astro-ph/9710029]; G. Gelmini and A. Kusenko, Phys. Rev. Lett. 82 (1999) 5202, [hep-ph/9902354]; G. Gelmini and 

A. Kusenko, Phys. Rev. Lett. 84 (2000) 1378, [hep-ph/9908276]; J.J. Blanco-Pillado, R.A. Vazquez and E. Zas, Phys. Rev. 

D61 (2000) 123003, [astro- ph/9902266]; J.L. Crooks, J.O. Dunn and P.H. Frampton, (2000), astro-ph/0002089. 
[4] for a review, see F. Halzen, in Boulder 1998, Neutrinos in physics and astrophysics, 524-569, [astro-ph/98 10368] and 

references therein. 

[5] J.G. Learned and S. Pakvasa, Astropart. Phys. 3 (1995) 267, [hep-ph/9405296]; F. Halzen and D. Saltzberg, Phys. Rev. 
Lett. 81 (1998) 4305, [hep-ph/9804354]; H. Athar, Astropart. Phys. 14 (2000) 217, [hep-ph/0004191]; S.I. Dutta, M.H. 
Reno and I. Sarcevic, Phys. Rev. D62 (2000) 123001, [hep-ph/0005310]; J. Alvarez-Muniz, F. Halzen and D.W. Hooper, 
Phys. Rev. D62 (2000) 093015, [astro-ph/0006027]; H. Athar, G. Parente and E. Zas, Phys. Rev. D62 (2000) 093010, 
[hep-ph/0006123]. 

[6] L. Bento, P. Keranen and J. Maalampi, Phys. Lett. B476 (2000) 205, [hep-ph/9912240]; H. Athar, M. Jezabek and O. 

Yasuda, (2000), [hep-ph/0005104]. 
[7] C. Lunardini and A. Y. Smirnov, Nucl. Phys. B583 (2000) 260 [hep-ph/0002152]. 

[8] H. Rang and G. Steigman, Nucl. Phys. B372 (1992) 494; GB. Larsen and J. Madsen, Phys. Rev. D52 (1995) 4282. 
[9] J. A. Adams and S. Sarkar, preprint OUTP-98-70P, and talk presented at the workshop The Physics of Relic Neutrinos, 
Trieste, September 1998; W.H. Kinney and A. Riotto, Phys. Rev. Lett. 83 (1999) 3366, [hep-ph/9903459] ; J. Lesgourgues 
and S. Pastor, Phys. Rev. D60 (1999) 103521 [hep-ph/9904411]. 
[10] P. de Bernardis et al, Nature 404 (2000) 955, [astro- ph/0004404] ; A.E. Lange et al, (2000), astro-ph/0005004; S. Hanany 

et al, (2000), astro-ph/0005123; A. Balbi et al, (2000), astro-ph/0005124. 
[11] J. Lesgourgues and M. Peloso, Phys. Rev. D62 (2000) 081301, (2000), [astro- ph/0004412]; 
[12] S. Hannestad, Phys. Rev. Lett. 85 (2000) 4203, [astro-ph/0005018]. 



31 



M. Orito, T. Kajino, G. J. Mathews and R. N. Boyd, astro-ph/0005446; S. Esposito et al., JHEP 09 (2000) 038, [astro- 
ph/0005571]; S. Esposito et al., (2000), astro-ph/0007419. 

M.J.S. Levine, Nuovo Cim. 48A (1967) 67; V.K. Cung and M. Yoshimura, Nuovo Cim. A29 (1975) 557; D.A. Dicus and 
W.W. Repko, Phys. Rev. D48 (1993) 5106, [hep-ph/9305284]. 

For recent discussions, see D.A. Dicus and W.W. Repko, Phys. Rev. Lett. 79 (1997) 569, [hep-ph/9703210]; A. Abbasabadi 
et al, Phys. Rev. D59 (1999) 013012, [hep-ph/9808211]; A. Abbasabadi, A. Devoto and W. W. Repko, Phys. Rev. D 63 
(2001) 093001 [hep-ph/0012257]. 

G.G. Raffelt, "Stars as Laboratories for Fundamental Physics: The astrophysics of neutrinos, axions, and other weakly 
interacting particles" , Chicago, USA: Univ. Pr. (1996) 664 V- 
See e.g. D. Notzold and G. Raffelt, Nucl. Phys. B 307 (1988) 924. 

S. Samuel, Phys. Rev. D 48 (1993) 1462; V. A. Kostelccky, J. Pantaleone and S. Samuel, Phys. Lett. B 315 (1993) 46; 

V. A. Kostelecky and S. Samuel, Phys. Rev. D 49 (1994) 1740, Phys. Lett. B 385 (1996) 159 [hep-ph/9610399]. 

For a discussion of the evolution equations see e.g. G. Sigl and G. Raffelt, Nucl. Phys. B 406 (1993) 423; and references 

therein. 

V. A. Kostelecky and S. Samuel, Phys. Rev. D 52 (1995) 621 [hep-ph/9506262]; S. Samuel, Phys. Rev. D 53 (1996) 5382 
[hep-ph/9604341]; J. Pantaleone, Phys. Rev. D 58 (1998) 073002. 
L. Stodolsky, Phys. Rev. D 36 (1987) 2273. 

G. Raffelt, G. Sigl and L. Stodolsky, Phys. Rev. Lett. 70 (1993) 2363 [hep-ph/9209276]. 

The interplay of collisions and refraction is illustrated also in A. Y. Smirnov, "New Aspects Of Neutrino Oscillations In 
Matter," In "Lea Arcs 1981, Proceedings, New and Exotic Phenomena" , 275-289. 

K. Enqvist, K. Kainulainen and M. Thomson, Nucl. Phys. B 373 (1992) 498; S. Dodelson and M. S. Turner, Phys. Rev. 
D 46 (1992) 3372; X. Shi, D. N. Schramm and B. D. Fields, Phys. Rev. D 48 (1993) 2563 [astro- ph/9307027]. 
J. Pantaleone, Phys. Lett. B 287 (1992) 128. 
E. Roulet, Phys. Rev. D47 (1993) 5247. 

See e.g. M.C. Gonzalez- Garcia et al., (2000), hep-ph/0007227. 
C. Lunardini and A. Yu. Smirnov, in preparation. 

M. Krumholz, S.E. Thorsett and F.A. Harrison, Ap. J. 506 (1998) L81, [astro-ph/9807117]; D.W. Hogg and A.S. Fruchter, 
Ap. J. 520 (1999) 54; IASSNS-AST-98-32. M. Schmidt, Ap. J. 523 (1999) L117, [astro- ph/9908206]; H.J.M. S. Mao, Astron. 
and Astrophys. 339 (1998) LI, [astro-ph/9908342]; see also E. E. Fenimore and E. Ramirez-Ruiz, (2000), astro-ph/9906125. 

E. Waxman and J. Bahcall, Phys. Rev. Lett. 78 (1997) 2292 [astro-ph/9701231]. 
T. Miyaji, G. Hasinger and M. Schmidt, MPE-478 [astro-ph/9910410]. 

F. W. Stecker and M. H. Salamon, Space Sci. Rev. 75, 341 (1996) [astro-ph/9501064]. 

P. Bhattacharjee and N.C. Rana, Phys. Lett. B246 (1990) 365; P. Bhattacharjee and G. Sigl, Phys. Rev. D51 (1995) 4079, 
[astro-ph/9412053]. 

U.F. Wichoski, J.H. MacGibbon and R.H. Brandenberger, (1998), hep-ph/9805419. 
P. Bhattacharjee, C.T. Hill and D.N. Schramm, Phys. Rev. Lett. 69 (1992) 567. 
See e.g. J. N. Bahcall, P. I. Krastev and A. Y. Smirnov, hep-ph/0103179. 



32 



